Mon. Not. R. Astron. Soc. 000, [THS] (2002) Printed 22 February 2012 



(MN I*TeX style file v2.2) 



(n: 



Probing the properties of Be star discs with 
spectroastrometry and NLTE radiative transfer modelling: 
CMi^ 



H.E. Wheelwright^'^t J.E. Bjorkman^, R.D. Oudmaijer^, A.C. Carciofi'^, 
X)" K.S. Bjorkman^ and J.M. Porter^J 

^ Max- Planck- Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 
P^H ■ ^School of Physics & Astronomy, University of Leeds, Woodhouse Lane, Leeds LS2 9JT, UK 

^ ■ ^Department of Physics & Astronomy, University of Toledo, MSlll 2801 W. Bancroft Street, Toledo, OH 43606, USA 

, "^Instituto de Astronomia, Geofisica e Ciencias Atmosfericas, Universidade de Sao Paulo, Rua do Matdo 1226, Cidade Universitdria, 
. 05508-900 Sao Paulo, Brazil 

l— I" ^Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Egerton Wharf, Birkenhead CH4I 1 LD, UK 

00 



Accepted yyyy Month dd. Received yyyy Month dd; in original form yyyy Month dd 



o 

^ ; ABSTRACT 

^ ■ While the presence of discs around classical Be stars is well established, their 

origin is still uncertain. To understand what processes result in the creation of these 
discs and how angular momentum is transported within them, their physical proper- 
ties must be constrained. This requires comparing high spatial and spectral resolution 
data with detailed radiative transfer modelling. We present a high spectral resolution, 
■ i?^80,000, sub milli-arcsecond precision, spectroastrometric study of the circumstellar 

disc around the Be star /3 CMi. The data are confronted with three-dimensional, NLTE 
. radiative transfer calculations to directly constrain the properties of the disc. Further- 

PsJ ■ more, we compare the data to disc models featuring two velocity laws; Keperian, the 

I prediction of the viscous disc model, and angular momentum conserving rotation. It 

04 ■ is shown that the observations of /3 CMi can only be reproduced using Keplerian ro- 

tation. The agreement between the model and the observed SED, polarisation and 
spectroastrometric signature of /3 CMi confirms that the discs around Be stars are 
well modelled as viscous decretion discs. 
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1 INTRODUCTION 

Be stars are non-supergiant stars with a spectral type B that 
either exhibit Hi Balmer emission lines in their optical spec- 
tra or have done so in the past ((Collins ,1987, ). It has been 
thought for some time that this Bal mer line emi ssion origi- 
nates in gaseous disc- like structures (|Struvell93ll ). The pres- 
ence of rotating, flattened material, i.e. a disc, was initially 
inferred from the typical double-peaked emission line pro- 
file that is exhibited by many Be stars. This has since been 
confirmed by long-baseline optical and near-infrared inter- 

* Based on observations conducted at the European Southern 
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082.D-0140. 

f E-mail: hwheelwright@mpifr-bonn.mpg.de 
I Deceased 



ferometry which has resolved small scale d iscs around several 
Be stars (see e.g. iQuirrenbach et al.|[l"997l ). However, despite 
the progress made in recent years, man y questions remain 
with regard to the origin of these discs l|Porter fc Riviniua 
l2003t ). 

Most Be stars rotate rapidly. Indeed, it appears that Be 
stars preferentially rotate at a significant fraction, ~80 per 
cent or more, of the critical ve locity where the centrifugal 
force and gravity are balanced (|Porteij|l996l : iRivinius et al.l 
12006). This rapid rotation may allow material to be lifted 
off the surface of the star and into orbi t, thus forming the 
gaseous circumstellar disc (|Struvelll93il ). However, this re- 
quires the star to rotate perilously close to its critical veloc- 
ity, and it is not clear whether Be stars rotate this r apidly or 
not (see e.g. iTownsend et al.|[20o3 : ICranmeill2005l ). There- 
fore, alternative mechanisms such as non-radial pulsations 
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Rivinius et alj 20031: [Cranmei 200£ ) and magnetic fields 
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2OO2I : iBrown et al 



20081 ) have been pro- 



posed as a viable mechanism to lift material ofi' the central 
star. 

Regardless of how material is ejected into orbit, another 
mechanism is required to distribute the material through- 
out the disc. T he viscous decretion disc model (VDDM, 
iLee et aD Il99ll ) features angluar momentum (AM) trans- 
port by turbulent viscosity to lift material into higher or- 
bits, thereby causing the disc to grow. A necessary re- 
sult of this model is that material is in Keplerian rotation 
throughout the disc. Other models, that do not have an 
AM transport mechanism, such as non-viscous magnetically 
confined discs and the origina l wind compressed disc model 
l|Biorkman fc Cassinellil Il993l ') , exhibit angular momentum 
conserving (AMC) kinematics. Here we directly confront 
these two fundamental predictions with astrometric data 
and physical modelling. 

Evidence that the discs around Be stars are well mod- 
elled as viscous decreti on discs is mounting . An example is 
provided by the work of lCarciofi" et all l|2009l ). who modelled 
the cyclic variabili ty of C Tau wit h perturbations in a viscous 
disc. The model of lCarciofi et al.! (2009) reproduces both the 
VLTI/ AMBER observations of ^ Tau and the temporal vari- 
ations of its Hq and Br7 emission and thus provides strong 
support for the viscous disc scenario. As another example, 
Keplerian rotation, a feature of the VDDM, is suggested 
by m any Be emission line profiles jHummel fc VranckenI 
I2OO0I ). As a result, the Keplerian viscous disc is currently 
the favoured model o f Be s tar discs. However, as noted by 
iHummel fc VranckenI (|2000l ), based on the line profiles alone, 
it is difficult to conclusively discount alternate velocity laws. 

A number of discs around Be stars have now 
been directl y probed with optical/ NIR inter ferometry (see 
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Vakih et al. 1998; Mcillan d et al.l l2 007a b: 

2009 ; Pott et al, ^2010; Delaa et all ^011; 

20 111 : iKraus et al.l l2012l : ?). Several of these studies also 
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find evidence for discs in Keplerian r otation. However, there 
are a few notab le exceptions (see iMeilland et al.l l2007al : 
iDelaa et al."201lV Therefore, we have yet to arrive at a com- 
plete understanding of the properties of Be star discs. Fur- 
thermore, many of the earlier studies that employed spectro- 
interferometric data were limited in terms of spectral reso- 
lution and/or had a sparse u,v coverage. Thes e limitations 
have been overcome in a few cases (see e.g. iKraus et al.l 
I2OI2I ). Nonetheless, few studies confront spectrally and spa- 
tially resolved data with detailed three-dimensional, NLTE 
radiative transfer calculations. Here we address this issue 
and confront unique high spectral resolution, sub milli- 
arcsecond (mas) precision spectroastrometric data with a 
detailed model of Be star discs. 

Spectroastrometry is a complementary technique to 
spectro-interferometry that offers a valuable insight into the 
kinematics of unresolved structures. The technique utilises 
the spatial information present in a long-slit spectrum to 
deli verer sp atial informati on with sub-mas precision (see 
e.g. lBaileyf ll998: Whclan fcGarcial I2OO8I. Spectroastrom- 
etry can offer a higher spectral resolution than spectro- 
interferometry and offers an efficient way to probe circum- 
stellar kinematics on small scales. Therefore, such data are 
well suited to comparison with model calculations. 



We have observed a sample of Be stars with spectroas- 
trometry. Here we present our results on /3 CMi, a B8 type 
star located at a distance of 52 pc. The disc around this 
object was recently resolved with the VLTI and CHARA 
interferometers as reported by iKraus et all l|2012l ). Using a 
kinematical model, these authors show that their observa- 
tions are best described assuming Keplerian rotation. Here, 
we fit our high spectral resolution data and photometric and 
spectro-polarimetric data with a state-of-the-art NLTE ra- 
diative transfer model to probe the physical properties of 
the disc. 



2 OBSERVATIONS AND DATA REDUCTION 

P CMi was observed on the 9th of Decem ber 2008 with the 
UVES spectrograph ijDekker et all I2OO0I ) mounted on the 
VLT-UT2. The red arm of the instrument was used as at- 
tention was focused upon the Ha line. Observations where 
conducted using the 600 g/mm grating and the Ha filter. 
The MIT-LL CCD was employed and the spatial pixel size 
of 0.16" ensured the seeing, which was ~0.8", was well sam- 
pled. Observations were conducted with a slit width of 0.5", 
which resulted in a spectral resolution ~80, 000 or 4 kms~^. 
The spectral range was ~ 6545-6580 A. Data were obtained 
at four different slit position angles, 0, 90, 180 and 270° . This 
is to identify, and eliminate, any syste matic artifacts in 
the^spectroastrometric signatures (see e.g. iBrannigan et al.] 
[200 61" ). Spectra of a ThAr lamp were used to wavelength cal- 
ibrate the data. 

Data reduction was conducted using IRAF and routines 
written in idl. Flat field frames, which were first corrected 
for the average bias level, were combined. The averaged fiat 
field was then normalised. Finally, raw data were corrected 
using an average bias frame and the normalised, average 
flat frame. The total intensity long-slit spectra were then 
extracted in a standard fashion. 

Spectroastrometry was performed by fitting Gaussian 
profiles to the spatial profile of the long-slit spectra at each 
pixel along the dispersion axis. This allowed the centroid 
position to be determined as a function of velocity. The con- 
tinuum position exhibited a trend in the dispersion direction 
which was removed by fitting a 1-D polynomial function to 
it. The average positional spectra for anti-parallel position 
angles were combined to form the average North-South (NS) 
and East- West (EW) positional spectra. All such spectra 
were assessed visually to exclude features only present at a 
single position a ngle (for details on the spectroastrometric 
techn ique see e.g. lOudmaiier et al.ll2008l : IWheelwright et ahl 
I2OIOI ). 

The spectroastrometric signature of /? CMi is shown in 
panel d) of Figure [T] As a result of the high SNR, the av- 
erage positional precision (the rms of the centroid in the 
continuum) is less than 0.2 mas. A clear spectroastrometric 
signature is o bserved and appears as expected for a rotating 
disc (see e.g. lOudmaiier et al.ll201ll ). The amplitude of the 
excursions in the EW and NS directions is of order 1 mas, 
and allows us to derive the Position Angle (PA) traced by 
the major axis of the disc of 130±5°. This is consistent with 
the polarisation angle of ~ 45°, as this is expected to perpen- 
dicular to the disc. The calculated PA is also ess entially con- 
sistent with the value of 140 ±1.7° reported bv lKraus et al.l 



Spectroastrometry and NLTE modelling of /3 CMi 's disc 3 



1I2012I ). We note that the spectroastrometric signature ex- 
hibits no sign of an additional component suc h as a polar 
outflo w. This is consistent with the finding of iKraus et alJ 
l|2012l) who report that the complex differential phase sig- 
natures of this star are, in-part, the result of the visibility 
function passing through a null. 



point source cata logue jCutri et al. 20o3 'l. t he IRAS poin t 
source catalogue (|Beichman et al.l 19881 ') and l|Ducatill2002l '). 
The UV spectrum was obt ained from the INES databasqj 
(see lWamsteker et al.lboool '). and broad-band measurements 
of the UV and optica l flux were taken from I Johnson et al.l 
(|l966l ): lOucatil (|2002l '). Optical spectra were taken from the 
HPOL data-base mentioned above. 



3 MODELLING 

The main aim of this work is to directly constrain the 
physical properties and kinematic structure of the disc. 
We do this by comparing the spatial and kinematic infor- 
mation contained within the spectroastrometric data and 
archival photometric and spectro-polarimetric data to the 
state-of-thc art radiative trans fer code HDUST developed 
by ICarciofi fc BiorkmanI l|2006l ). This is done using using 
two rotation laws: Keplerian rotation, v{r) cx r~2, and the 
AMC case, v{r) oc r~^. 



3.1 HDUST: NLTE radiative transfer 

HDUST is a fully 3D, non-local thermodynamic equilibrium 
(NLTE), Monte Carlo code which solves the problems of ra- 
diative transfer, radiative equilibrium and statistical equilib- 
rium for pre-defined gas density and velocity distributions. 
Details of t he code can be found in ICarciofi fc BiorkmanI 
l|2006l . l200d ). Here, the code is used to generate the SED, 
polarisation and Ha line profile of a rotating star surrounded 
by a gaseous equatorial disc, in addition to making images 
of the system. We restrict ourselves to the application of the 
code and refer the reader to the references above for more 
details. 

Determining the self consistent structure of Be star discs 
is a complex problem, and requires solving the equations 
for vertical hydrostatic equilibrium an d radial viscous diffu- 
sion in non-isothermal co nditions (see Carciofl fc BiorkmanI 
l2008l 'l. This was done bv lCarciofl fc BiorkmanI (120081 ) who 
demonstrated that in a low density viscous disc, the temper- 
ature structure is approximately isothermal vertically. Such 
a disc is well described by a model in which the vertical den- 
sity follows a Gaussian distribution and the density exhibits 
power-law depencence in the radial direction. (3 CMi exhibits 
a low polarisation, indicative of a low density disc. Conse- 
quently, we use this analytical description for the structure 
of the disc. 



3.2 Archival data 

Spectro-polarimetric data were obtained using the HPOL 
spectropolarimeter, mounted on the Pine Bluff Observatory 
telescope. j3 CMi was observed on the night of 22/04/1991 
using a dual Reticon array detector spanning the wavelength 
range of 3200-7600 A and with a spectral resolution of 25 A. 
The Reticon detector was replaced with a CCD detector and 
two new gratings which extended the wavelength coverage to 
3400- 10500 A and improved the spectral resolution to ~ lOA. 
P CMi was then observed several times with HPOL (1995 fc 
2000). Comparisons of different observations indicated that 
the flux blue-wards of the Balmer jump is possibly affected 
by systematics. Infrared fluxes were taken from the 2MASS 



4 RESULTS 

4.1 Modelling procedure and results 

The initial step in reproducing the observations was deter- 
mining the properties of the central star. We estimated the 
radius of the star and the reddening towards it by reproduc- 
ing the observed UV and optical SED with an atmospheric 
model appropriate for the spectral type of P C Mi (B8Ve, 
log = 4.0, Teff = 12, 000 A'l. lFremat et al.l(|2005l ) report val- 
ues of vsini = 231 kms~^ and Vcrit = 380 kms~^. Allowing 
for a ratio «rot/wcrit = 0.8 — 0.85 suggests the inclination is 
i ^ 45 — 50°. in app roxi mate agreem e nt wi th the results of 
iTvcner et aU (|2005l ) and lKraus et all (|2012l '). We assume an 
inclination of i = 40 — 45° . 

It was assumed that the star rotates relatively rapidly, 
at approximately 80-85 per cent of its critical velocity. Grav- 
ity darkening effects due to the rapid rotation are taken into 
account in HDU ST using the von Zeipel flux distribution 
llvon Zeipej|l924l '). The ratio between the equatorial and po- 
lar radii and temperature were determined based on the stel- 
lar rotation and critical velocity. To conduct the modelling, 
we explored the available parameter space deflned by the 
ranges above. This was done varying i, Kot, Vcrit and i?niax 
by hand as the time required to iteratively fit all the obser- 
vations was prohibitively large. Vxurb and po were treated as 
free parameters. We note that reproducing the photometric, 
polarimetric and spectroscopic data rules out many parame- 
ter combinations that might represent local minima if fewer 
types of data were considered. 

For the Keplerian model, we adopted an inclination of 
i — 40° and rotation at 85 per cent of the critical velocity 
as this allows a better match to the observed line profile. 
This does not produce the best fit to the HPOL SED, but 
difference between the predicted SED and that observed is 
within the uncertainty in the flux measurements. The in- 
frared excess of jS CMi is comparatively small. Reproducing 
this requires a low disc density, which is consistent with the 
low polarisation of the star. The best fltting density was 
deflned as the value that reproduced the line-to-continuum 
ratio of the Ha emission. In principle, the central reversal 
can be recreated by increasing the line-width. However, this 
increases the line-to-continuum ratio. The density cannot be 
reduced to compensate without degrading the flt to the IR 
excess and thus the central reversal cannot be recreated ex- 
actly. We also note that the high velocity wings of the Ha 
emission cannot be reproduced. This is attributed to broad- 
ening that is not present in the model. 

For the AMC model, we adopted a rotation at 80 per 
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Table 1. The list of key parameters. 



Parameter 



ValueKcp ValucAMC Notes 



d (pc) 


52 


52 


1 


i° 


40 


45 


2 


Vrot/Vcrit 


0.85 


0.8 


3 


Rc (R©) 


4.68 


4.68 


4 




0.05 


0.05 


4 


Tp{ K) 


15,000 


15,000 


4 




25 


3.3 


5 


V-Tuib (kms-1) 


0.1 


14.5 


5 


M. (Mq) 


3.08 


3.97 


6 


Re/ Rp 


1.32 


1.27 


6 


Tp/Te 


1.48 


1.37 


6 


po (xlO~^^gcm~^) 


2.9 


4.6 


7 


L. (Lo) 


220 


220 


8 


i;sini (kms ^) 


182 


209 


9 


Vcrit(kms-i) 


333 


370 


10 



Notes: 1: Hippa rco s parallax llPerrvman et al.l Il997t ). 2: 
iTvcner et al] l l2005h fc iKraus et~al] ||2012|) . 3: fit SEP and line 
profile, 4: fit SED, 5: fit line profile, 6: Based on Rp and Vcrit, 
7: fit line flux and polarisation, 8: set, 9: c.f. 230 ili 16 km s~^ 
iFremat et al] ||2005|) and 210 ifc ~ 20 kms-l [Ab"t et al.l ||2002| '), 10: 
see e.g. iFremat et al] ||2005|) . 



cent of the critical velocity and an inclination of i = 45° 
as this allowed the Ha line profile to be reproduced. Since 
the velocity in the disc decreases more rapidly with radius 
than in the Keplerian case, the outer radius of the disc has 
to be significantly reduced to recreate the observed double- 
peaked line profile. To match the observed line-to-continuum 
ratio, the line-width also had to be increased. In this case, 
this resulted in a good fit to the central reversal. The fi- 
nal model line profiles, polarisation signatures, SEDs and 
spectroastrometric signatures are presented in Figure [T] and 
the associated parameters are listed in Table [T] It is clear 
that both scenarios are consistent with the entire polarimet- 
ric, photometric (including the IR photometry which is not 
shown) and spectroscopic data-set, preventing differentia- 
tion between the two velocity laws on this basis alone. 



4.2 Spectroastrometric results 

The crucial constraint on the disc properties is provided by 
the spectroastrometric data. The spectroastrometric signa- 
tures of the best fitting Keplerian and AMC models are pre- 
sented in panel d) of Figure [1] The amplitude of the spec- 
troastrometric signature of the AMC disc is smaller than 
that of the Keplerian model; due to the smaller size of the 
AMC disc. On the contrary, the signature of the Keplerian 
model reproduces the observed signature well (x^ ~1.5 ver- 
sus ~5). We emphasise that the model astrometric signa- 
tures were predicted by the models that best reproduce the 
observed SED, polarisation and line profile. The spectroas- 
trometric signatures were not fit in any way. Since only the 
Keplerian model recreates the observed spectroastrometric 
signature, the AMC scenario can be discarded. Therefore, 
we confirm that the kinematical structure of (3 CMi's disc is 
consistent with Keplerian rotation. 
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Figure 1. The final HDUST SEDs (a), optical polarisation (b) 
and line profiles (c) and the predicted spectroastrometric signa- 
tures (d). The observations (black) are presented alongside the 
final Keplerian (red) and AMC (blue) models. 



5 DISCUSSION AND CONCLUSION 

We present a comparison between photometric, polarimetric 
and spectroastrometric observations of the Be star P CMi 
and HDUST, a NLTE radiative transfer code developed to 
model Be stars and discs. Currently, the favoured model 
of the discs around Be stars is the viscous decretion disc 
model. We confront this model with spectrally resolved sub- 
mas precision observations which yield direct evidence of the 
kinematic structure of the disc around P CMi and provide 
a unique test of the viscous disc scenario. 
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We show that both Keplerian (the prediction of the 
VDDM) and AMC scenarios can reproduce the spectro- 
scopic, photometric and polarimetric data-set. However, we 
were able to subject the disc kinematics to a critical test 
via spectroastrometry which yields velocity resolved spatial 
information with sub-mas precision. Due to the rapid reduc- 
tion in rotational velocity with radius, the AMC disc that 
recreates the observed line profile must be smaller than the 
Keplerian disc. As a result, the spectroastrometric signa- 
ture of the best fitting AMC disc is smaller than that of 
the Keplerian model. Consequently, only the Keplerian disc 
recreates the spectroastrometric signature observed and the 
AMC scenario can be discounted. 

That the disc around /3 CMi can be modelled as a vis- 
cous disc in Keplerian rotation indicates that the transport 
of angular momentum within it is governed by turbulent vis- 
cosity. This addresses one of the two key questions regard- 
ing the discs of Be stars: how is material injected into the 
disc and how is angular momentum transported within it? 
Although we address the latter question, this does not con- 
strain the disc formation mechanism. It has been suggested 
that the disc around /3 CMi i s related to the n on-radial pul- 
sations this object exhibits l|Saio et al.l |2007| ). These data 
cannot directly assess this possibility. However, our mod- 
elling does yield the physical properties of the disc which 
may provide useful constraints to future models of disc for- 
mation. 

We conclude by emphasising that the comparison be- 
tween our extensive data-set and NLTE radiative transfer 
modelling provides a stringent test of the properties and 
kinematics of the disc of /3 CMi. The agreement between 
the model and observations in both the spectral and spatial 
domains clearly demonstrates that the discs around Be stars 
can be well modelled as viscous decretion discs. 
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